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Abstract. The question of obtaining well-defined criteria for multiple cri- 
teria decision making problems is well-known. One of the approaches dealing 
with this question is the concept of nonessential objective function. A cer- 
tain objective function is called nonessential if the set of efficient solutions 
is the same both with or without that objective function. In this paper we 
put together two methods for determining nonessential objective functions. 
A computational implementation is done using a computer algebra system. 

Key Words. Multiobjective optimization problems. Efficient (Pareto op- 
timal) solutions. Essential/nonessential objective functions. 

1 Introduction 

Multiple criteria decision making (MCDM) arise in connection with the so- 
lution of problems in the areas of economy, environment, business and engi- 
neering (see e.g. Refs. [UH]). Multiobjective programming is concerned with 
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the generation of solution sets for multiobjective problems that usually in- 
clude a large or an infinite number of points, known as efficient solutions. 
Those efficient points are then the candidates for an optimal solution of the 
MCDM problem. 

The question of obtaining well-defined criteria for MCDM problems is 
well-known. Often mathematical models done by inexperienced practitioners 
lead to redundant formulations, which are not only deceptive but also compu- 
tationally cumbersome (see Refs. 013 and therein). Sometimes the decision- 
maker end up without any decision support, while a simple reformulation of 
the problem would achieve the desired result. One of the approaches dealing 
with the question of obtaining well-defined criteria for MCDM is based on 
the concept of nonessential objective function. A certain objective function 
is called nonessential if the set of efficient solutions is the same both with or 
without that objective function. Information about nonessential objectives 
helps a decision maker to get insights and better understand a problem, and 
this might be a good starting point for further investigations or revision of 
the model. Dropping out nonessential functions leads to a problem with a 
smaller number of objectives, which can be solved more easily. For this rea- 
son, the identification of nonessential objectives is an important feature in 
the analysis of multiple criteria programs. 

The seminal papers by Gal and Leberling (Refs. El [9j) define and investi- 
gate nonessential objectives in linear multiobjective optimization problems; 
Gal and Hanne (Refs. El El) study the consequences of dropping nonessen- 
tial objectives functions in the application of MCDM methods. Recently, 
the concept of nonessential objective has been generalized by Malinowska to 
convex multiobjective optimization problems, and new definitions of weakly 
nonessential and properly nonessential objective functions were introduced 
and investigated (Refs. [TU | [TT|) : a new method to determine if a given objec- 
tive function of a certain linear problem is essential or not has been proved 



Here we put together, in a constructive and algorithmic way, the two 
available methods (Refs. El 1T21) for determining nonessential objective func- 
tions. A computational implementation is done using the computer algebra 
system Maple. The plan of the paper is as follows. Section [2] gives the nec- 
essary preliminaries and provides the notation used in the text. In Section [3] 
we develop the theory of nonessential objectives. Main result of the paper 
appears in Section HJ the algorithm to determine if a given objective function 
of a linear multiobjective problem is essential or not. Finally, in Section[S]we 
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provide some examples that show the applicability of our methodology and 
the convenience of the developed computer software. The paper ends with 
some conclusions, the references, and an appendix with all the Maple code 
that implements the proposed algorithm. 

2 Preliminaries and Notation 

In this section we present some general concepts and notations. We use 
superscripts for vectors (for example x 1 , or simply x when no confusion can 
arise), and subscripts for components of vectors (for example X\). All the 
vectors are assumed to be column vectors. The symbol 1 stands for the vector 
[1, . . . , 1] T . For two vectors x,x' EM. k we define the relations (Ref. [TBI) : 

x ^ x' Vi G {1, . . . , k} : Xi > x' i , 

x > x' O Vi G {1, . . . , : Xi > x\ A 3i G {1, . . . , k} : x { > x\ , 

x > x' 4^ Vi G {1, . . . , k} : Xi > x[ . 

Throughout this paper we study the following linear multiobjective opti- 
mization problem: 

max{F n+1 (a;) : x G X}, (1) 

where 

X = {x G R k : Ax ^ b , x ^ 0} , A G R mxk , b G M. m 
is the feasible set, and 

F n+1 (x) = Cx=[(c 1 ) T x,...,(c n+1 ) T xf , e eR fe (i = l,...,n+l), n>l 

is the vector of objective functions: fi : M. k —>■ R (i = 1, . . . , n + 1). We are 
using "max" to mean that we want to maximize all the objective functions 
simultaneously. This involves no loss of generality. In general it does not 
exist a solution that is optimal with respect to every objective function, and 
one is lead to the concept of Pareto optimality. 

Definition 2.1 A vector x° G X is said to be an efficient (Pareto opti- 
mal) solution of problem (T7]j if there exists no x G X such that F n+1 (x) > 
F n+1 (x°). The set of efficient solutions of problem (QP is denoted by X^ +1 . 



3 



3 Nonessential Objectives 



Let X% denote the set of efficient solutions of problem ([I]) without one ob- 
jective function k G {1, . . . , n + 1}. Without loss of generality we assume 
k — n+1. 

Definition 3.1 The objective function f n+ i is said to be nonessential in 
problem (TJP if X 1 ^ = X^ +1 . An objective function which is not nonessen- 
tial is called essential. 

We now recall two theorems that characterize a nonessential objective func- 
tion, and which will be used in the proposed algorithm to determine if a 
given objective function is essential or not. 

Theorem 3.1 (Ref. [2P The objective function f n+1 is nonessential in (TJP if 
the following holds: 

n 

c n+1 = ^a iC \ Oi>0 (i= l,...,n). (2) 

i=l 

Theorem 13.11 is very useful because condition (T5]) is easy to check. Unfortu- 
nately, it is only a sufficient condition. Theorem 13.21 below gives a necessary 
and sufficient condition for an objective function to be nonessential. 

Let X n+ i denote the set of solutions of the single objective optimization 
problem 

max{/ n+ i(a;) : x G X} . (3) 

In other words, 

X n+1 = {x° G X : Vx G X f n+1 (x°) > f n+1 (x)}. (4) 

Theorem 3.2 (Refs. [TU[ [77]) If the set X is nonempty and bounded (in 
other words, if X is a convex polyhedron) , then the objective function f n+ x is 
nonessential in (CP if and only if the following three conditions hold: 

(i) Vx G X\Xl 3x' G R k : F n+l (x') > F n+1 (x); 

(ii) XI n X n+1 ^ 0; 

(Hi) X%CX« +1 . 
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4 Main Results 

The theory described in the previous section enable us to work out on a 
computational algorithm to test if an objective function f n+ i in problem 
(CD) is essential or not. The proposed algorithm consists of eight steps and 
each one has been implemented in the computer algebra system Maple (see 
Appendix on page [19]) . 

Step 0. Does c n+1 = £" =1 > (i = 1, . . . , n)? 

If the answer is "TRUE", then the objective function / n+1 is nonessential by 
Theorem 13. 11 Otherwise, we go to Step 1. 

We implement Step as a Maple command glm (Gal-Leberling method), 
which receives in its first argument a list with the objective functions, and 
in its second argument the number of variables of the problem. 

Example 4.1 Consider the multiobjective optimization problem 

max{F 4 (x) : x E X}, 

where 

fi(x) = xi + 3x 2 , 
f 2 (x) = 3xi, 
f 3 (x) = 2xi + x 2 , 
f 4 (x) = -3xi - x 2 , 

and IcK 2 . With our Maple package we do: 

> glm( [xl+3*x2,3*xl,2*xl+x2,-3*xl-x2] ,2) ; 
false 

We conclude that the answer to Step is "FALSE" , so we go to Step 1. 
Now, let us change the order of objective functions as follows: 

fi(x) = xi + 3x 2 , 
f 2 (x) = 3xi, 
/ 3 (x) = -3xi - x 2 , 
/ 4 (x) = 2xi + x 2 . 

This time our g Im procedure gives 
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> glm( [xl+3*x2,3*xl,-3*xl-x2,2*xl+x2] ,2) ; 



true 



meaning that the answer to Step is "TRUE". Thus, one can conclude that 
the objective function fi{x) = 1x\ + x 2 is nonessential. 

Example 4.2 Consider the problem: 



and IcK 3 . We have: 

> glm( [xl+x2+x3,-xl+x2+x3,xl+x2] ,3) ; 
false 

Answer to Step is "FALSE" , and we go to Step 1. 

Step 1. We test condition (i) of Theorem 13.21 

Our method is based on the following observations. Let 

U = {x G R k : Cx > 0}. 

Remark 4.1 (Ref. IfU^fy, then condition (i) of Theorem[3lM holds. 

Theorem 4.1 (Ref. \TQ) A sufficient condition for X^ +1 = X is U = 0. // 
intX ^ (where int stands for the interior of a set), then this condition is 
also necessary. 

In Step 1 we solve the problem: does U 0? If the answer is "FALSE", then 
Xg +1 = X and we go to Step 2. Otherwise, we know that condition (i) of 
Theorem 13.21 holds and we go to Step 5. In order to verify equality of sets 
U = we solve the problem: 



max{F 3 (x) : x E X}, 



where 



h{x) 



Xi + x 2 + x 3 , 
—Xi + x 2 + x 3 
X\ + x 2 , 





where 



V = {(x,v) e E' 



k+n+l . 



Cx + v = , v ^ 0} . 
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Remark 4.2 (Ref. [X^ ) One has U = if and only if problem /ias zero 
as the optimal objective function value. 



Example 4.3 Consider the problem from Example \4-1 

max{F 4 (x) : x E X}, 

where 

/i(x) = Xi + 3x 2 , 
/ 2 (a?) = 3a;i, 
/ 3 (x) = 2xi + x 2 , 
/4(x) = -3x x - x 2 , 

and X C M 2 . Using our Maple command stepl we obtain: 

> stepl (2, [xl+3*x2, 3*xl, 2*xl+x2, -3*xl-x2] ) ; 
false 

We conclude that Step 1 has answer "FALSE" . Thus, X% = X and we go 
to Step 2. 



Example 4.4 Consider again the problem from Example 4-2: 

max{F 3 (x) : x G X}, 

where 

f x {x) = Xi +x 2 + x 3 , 
h{x) = -xi + x 2 + x 3 , 
/ 3 (x) = Xi +x 2 , 

and IcK 3 . We obtain 

> stepl (3, [xl+x2+x3, -xl+x2+x3, xl+x2] ) ; 
true 

so the answer to Step 1 is "TRUE". Thus, condition (i) of Theorem \3.2\ 
holds and we go to Step 5. 
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Step 2. Let U = {x G R k : Cx > 0}, where C = [(c 1 ) T , . . . , (c n ) T ] T . In 
Step 2 we address the following question: does [7^0? 

The method we use is the same as the one described in Step 1. If the 
answer is "FALSE", then Theorem 14.11 implies = X and the objective 
function f n+ i is nonessential. Otherwise, we go to Step 3. 

Example 4.5 Let us consider again the problem from Example \4-3\ Using 
our Maple command step2 one has: 

> step2(2, [xl+3*x2, 3*xl, 2*xl+x2, -3*xl-x2] ) ; 

true 

Since the answer to Step 2 is "TRUE" we go to Step 3. 

Step 3. In this step we solve the problem: does intX ^ 0? 
Our method is based on the following remark. 

Remark 4.3 If intX ^ 0, then problem max{a : (x,v,a) G V}, where 

V = {{x,v,a) E R k+m+1 : Ax + v + al = b,v^0,x^el,a>0} 

with e > sufficiently small, has an optimal objective function value greater 
than zero. 

If in Step 3 the answer is "TRUE", then Theorem 14. II implies X% ^ X, and 
the objective function f n+ i is essential. Otherwise, we go to Step 4. 

In order to use the simplex package already available from the Maple 
system, we put e = 0, 001. We note that by default we are assuming that all 
x variables are greater or equal than zero (the user does not need to mention 
this explicitly in the definition of the set X while using our Maple package). 

Example 4.6 Let us continue the problem from Examples \4-3\ and \4-5\ with 

X = {x G R 2 : xi < 1, x 2 < 1, x ^ 0}. 

Using our Maple command step3 we obtain 

> step3(3,{xl <= 1, x2 <= 1}); 

true 
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Since the answer to Step 3 is "TRUE", the objective function f^{x) = — 3xi~ 
x 2 is essential. 

Now, let us consider a different problem by changing the set X as follows: 



Since Step 3 has the answer "FALSE", we go to Step 4- 

Step 4. In this step we solve the problem: does X 1 ^ = XI 

Our method is the following: we compute a vertex x° of X and test if x 1 
is an element of X^. 

Theorem 4.2 (Ref. [73]) Let x° G X be given. Solve the problem 



S = {(x,e)E R k+n :xeX, f t (x) - e t = /.(x ), e 4 > 0, i = 1, . . . , n} . 



The vector x° is efficient if and only if problem has zero as the optimal 
objective function value. 

If x° is not efficient, then the answer from Step 4 is "FALSE" and the 
objective function / n+1 is essential. Otherwise, we compute the next vertex 
of X and check if it is efficient or not. Our procedure stops as soon as a 
non-efficient vertex is found. 

Example 4.7 Consider the problem: 



X = {x e M 2 : Xi + x 2 < 1, -Xi - x 2 < ^ 0}. 



Now we obtain 



> step3 (2 , {xl+x2<=l , -xl-x2<=-l}) ; 
false 





with 



max{F 3 (x) : x E X}, 
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where 

/ 2 (x) = X 1? 

/ 3 (x) = -3x 1 - x 2 , 

and 

X = {x eM 2 : Xi+x 2 <1, -Xi - x 2 < -1, x ^ 0}. 

Answers from Steps 1 to 4 ar & easily obtained from the respective commands 
of our Maple package: 

> stepl(2, [xl + x2, xl, -3 xl - x2] ) ; 

false 

> step2(2, [xl + x2, xl, -3 xl - x2] ) ; 

true 

> Step3(2,{xl+x2<=l,-xl-x2<=-l}) ; 

false 

> step4( [xl+x2,xl,-3*xl-x2] ,{xl+x2<=l,-xl-x2<=-l},2) ; 

false 

We conclude that fs(x) = —3x\ — x 2 is essential and that X\ C X%. 

In the case all vertices are efficient, two situations may appear: objective 
function may be essential (answer "FALSE") or not (answer "TRUE"). 
To distinguish between these two cases we apply the following remark. 

Remark 4.4 Let x 1 , x 2 , . . . , x p be all vertexes of X and intX = 0. Then, 
Xg = X if and only if there exists a vector w > with Y^7=i Wi = 1 suc h 
thatw T F n (x 1 )=w T F n (x 2 ) = --- = w T F n (x p ). 

Proof. As far as intX = 0, we have X C H, where if is a hyperplane. 
Therefore, X^ = X if and only if there exists w > (X)iLi w % = 1) such that 
for all x in X, x is a solution of the problem m&x{w T F n (x) : x G X} (see for 
instance Ref. [EH p. 54). This completes the proof. □ 

In order to use the simplex method as provided by Maple, we change 
condition from Remark 14.41 into the form 

n 

3w E W 1 : w T F n (x 1 ) = w T F n (x 2 ) = ■■■ = w T F n (x p )^Wi = l,w T l > el, 
where s > is sufficiently small. In our procedure we set e = 0, 00001. 
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Example 4.8 Let us continue the problem from Examples \4-S\ \4-5\ and \4-6] 

with 

X = {x G M 2 : xi + x 2 < 1, -Xi - x 2 < -1, x ^ 0}. 
Our Maple command step4 give us 

> step4( [xl+3*x2,2*xl+x2,3*xl,-3*xl-x2] ,{xl+x2<=l ,-xl-x2<=-l},2) ; 
true 

27ms, X\ = X and it follows that f^{x) = —3xi — x 2 is nonessential. 

Now we show an example where all vertices are efficient but X^ ^ X. 
Example 4.9 Let us consider the problem 

max{F 3 (x) : x G X}, 

where 

fi{x u x 2 , x 3 ) = -x\ - 2x 2 + 2x 3 , 
f 2 (x 1 ,x 2 ,x 3 ) = 2x 1 + 3x 2l 
/ 3 (xi, x 2 , x 3 ) = -xx -x 2 - 2x 3 , 

and 

X = (x G IR 3 : x 2 + x 3 < 2, —x 2 — x 3 < —2, x\ + x 2 + x 3 < 3, 

— xi — x 2 — x 3 < —2, x\ + x 2 < 2 , x ^ o| . 

Using our Maple commands we obtain: 

> stepK [-xl-2*x2+2*x3, 2*xl+3*x2,-xl-x2-2*x3] ) ; 
false 

> step2 ( [-xl-2*x2+2*x3 , 2*xl+3*x2 , -xl-x2-2*x3] ) ; 
true 

> step3 (3 , {x2+x3<=2 , -x2-x3<=-2 , xl+x2+x3<=3 , -xl-x2-x3<=-2 , xl+x2<=2}) ; 
false 

> step4( [-xl-2*x2+2*x3,2*xl+3*x2,-xl-x2-2*x3] , 

{x2+x3<=2 , -x2-x3<=-2 , xl+x2+x3<=3 , -xl-x2-x3<=-2 , xl+x2<=2} , 3) ; 

false 
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We conclude that the objective function f% is essential and that X 2 E C X\. 



Step 5. In this step we calculate all vertexes of X n+ i (see @). 

Let X^ +1 = {x 1 , x 2 , . . . , x q } be the set of all vertexes of We have: 

x EM. k : x = ctjxi , otj — 1 , ay > 0, j — 1, . . . , q 

i=i i=i 



Example 4.10 Consider again the problem from Example \4-4\ with 

X = {x eR 3 : X! <l,x 2 <l,x 3 <l,x^ 0}. 
Our procedure step5 give us 

> step5( [xl+x2+x3,-xl+x2+x3,xl+x2] , {xl<=l , x2<=l , x3<=l}) ; 

{{xj = 1, x 2 = 1, x s = 1} , {xt =l,x 2 = 1, % = 0}} 

Hence, Xf = {[1, 1, 0] T , [1, 1, 1] T } and 

X 3 = {x:x = a[l,l,0] T + (1 - a)[l, 1, 1] T ,0 < a < 1}. 

Step 6. In this step we solve the following problem (condition (ii) of 
Theorem does X n+l nl^§? 

The basic idea of our method is to use Theorem 14.31 
Theorem 4.3 (Ref Let Z = {x\ x 2 , . . . , x q } C X. If Z C X\X%, then 

q f/ "| 

xGK t :i = ^ ctjX\ y~] aij ; = 1, otj > 0, j = 1, . . . , q > C X \ . 
j=i i=i J 

Having in mind Theorem l4.3t it is sufficient to consider only vertexes of X n+1 . 
Applying Theorem 14. 21 we check if there exists a vertex of X n+ i which belongs 
to X E . If the answer is "TRUE" , then the condition (ii) of Theorem 13.21 holds 
and we go to Step 7. Otherwise, we conclude that the objective function f n+ i 
is essential. 



Example 4.11 Let us continue the problem from Examples 4-4 an d 4-10 
Our procedure step6 give the answer "TRUE" 
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> step6 ( [xl+x2+x3 , -xl+x2+x3 , xl+x2] , {xl<=l , x2<=l , x3<=l}) ; 

true 

and we proceed to Step 7. 

Example 4.12 Now we consider a problem borrowed from Ref. \4\- 

max{F 3 (x) : x G X}, 

where 

fi(x) = Xi + 2x2 — X 3 + 3X4 + 2^5 + £7, 
f 2 (x) = x 2 + x 3 + 2x 4 + 3x 5 + x 6 , 
f 3 (x) = Xx + x 3 - X4 - x 6 - x 7 , 

and 

X = {x E R 7 : xi + 2x 2 + x 3 + x 4 + 2x 5 + x 6 + 2x 7 < 16, 
— 2xi — x 2 + X4 + 2x 5 + x 7 < 16, 
— x\ + x 3 + 2x5 — 2x 7 < 16, x 2 + 2x 3 — X4 + £5 — 2x6 — X7 < 1, x ^: 0}. 

With our Maple command step6 we obtain 

> step6 ( [xl+2*x2-x3+3*x4+2*x5+x7 , x2+x3+2*x4+3*x5+x6 ,xl+x3-x4-x6- 

{xl+2*x2+x3+x4+2*x5+x6+2*x7<=16 , -2*xl-x2+x4+2*x5+x7<=16 , 
-xl+x3+2*x5-2*x7<=16 , x2+2*x3-x4+x5-2*x6-x7<=l}) ; 

false 

and since the answer is "FALSE", we conclude that the objective function 
f 3 (x) = Xi + x 3 — x 4 — x 6 — x 7 is essential. 

Step 7. In this step we solve the following problem (condition (iii) of 
Theorem [32D: does X E C 

Our method is based on the following observations. 

Proposition 4.1 (Refs. QJJ, IT8\) If the vector-valued function F n is one-to- 
one on the set X E; then condition (iii) of Theorem \3.2\ holds. 

Let 

(X n E ) = {x i -x> :x\x ] EX n E }. 
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Lemma 4.1 (Ref. 19) The vector-valued function F n is one-to-one on the 
set X E if and only if KerF n H (X E ) = (Ker stands for the kernel of a 
map). 

In practice it is usually impossible to determine the set (X E ). For this reason, 
in our Maple procedure we use the following set: 

where x l ,x° are vertexes of X E and x° is free but fixed. Obviously, (X E ) C 

Remark 4.5 If KerF n D Lin{(X^ VE )} = ; £/ien condition (Hi) of Theo- 
rem\3M holds (that is, X E C 

If the answer from Step 7 is "TRUE", then the objective function is 
nonessential. Otherwise, we know that X E +1 C X]^. 



Example 4.13 Let us continue the problem from Examples \4-4i \4-10{ an d 
\4-H\ We obtain: 

> step7 ( [xl+x2+x3 , -xl+x2+x3 , xl+x2] , {xl<=l , x2<=l , x3<=l} , 3) ; 
true 

The answer from Step 7 is "TRUE", hence the objective function fs(x) = 
X\ + x 2 is nonessential. 



5 Illustrative Examples 

We have implemented all the steps described in Section 4 together in a single 
Maple command called nonessential (see Figured]). This main procedure 
receives a list with the definition of the objective functions, and the set 
X of constraints in the second argument. Below we give some examples 
of computer sessions with our Maple package. The interested reader may 
download it from [http://www.mat.ua.pt/delfim/essential7html]] and 
find there more examples than the ones we are able to provide here. We invite 
and welcome the reader to experiment our Maple package with her/his own 
problems. 
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Step 



fn+l I s 

nonessential 



Step 1 



Step 2 



fn+l IS 

nonessential 



fn+i is 
essential 



fn+i is 
essential and 

E i= A £ 



Step 3 



Step 4 



/n+l is 

nonessential 



fn+l is 

essential 



fn+i is 
nonessential 



Step 5 



Step 6 



Step 7 



X™ +1 C X n E 



Figure 1: Scheme of the method implemented in Maple. 

We begin with a simple example with three objective functions and three 
variables. 

> nonessential ( [xl+x2,xl+x2+x3,-3*xl-3*x2-x3] ,{xl+x2+x3 <= 1}); 
Objective function -3*xl-3*x2-x3 is essential from Step 3 

Next we consider a problem with four objective functions and five vari- 
ables. It turns out that the problem may be reduced to a simpler one with 
the same set of efficient solutions. 

> st := timeO : 

nonessential ( [xl+x2+x3+x4+x5 , -xl+x2+x3+x4+x5 , -xl-x2+x3+x4+x5 , xl+x2] , 

{xl<=l ,x2<=l , x3<=l , x4<=l ,x5<=l}) ; 
printf("%a seconds\n" ,time () - st) ; 



Objective function xl+x2 is nonessential from step 7 
2.072 seconds 
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The only situation where our Maple procedure nonessential can not 
conclude that a given function is essential or not, is when one reaches Step 7 
and the answer is not true. 

> st := timeO : 

nonessential ( [-x3-x4,-x5-x6,-x4-x6] , 

{xl+3*x2<=24 , 3*xl+x2<=24 , xl+4*x2+x3-x4<=40 , 
-xl+4*x2-x3+x4<=-40 , 4*xl+x2+x5-x6<=40 , -4*xl-x2-x5+x6<=-40}) ; 
printf("°/„a seconds\n" ,time () - st) ; 

X_E~3 C X_E~2 from step 7 
7.251 seconds 

Our Maple package is useful to identify redundant objective functions. 
We finish with an example where the mathematical model can be simplified 
by elimination of two of the objective functions. 

> nonessentiaK [xl+3*x2,2*xl+x2,3*xl,-3*xl-x2] , 

{xl+x2<=l,-xl-x2<=-l}) ; 
Objective function -3*xl-x2 is nonessential from step 4 

> nonessential ( [xl+3*x2 , 2*xl+x2 , 3*xl] , {xl+x2<=l , -xl-x2<=-l}) ; 
Objective function 3*xl is nonessential from step 7 

> nonessential ( [xl+3*x2 , 2*xl+x2] , {xl+x2<=l , -xl-x2<=-l}) ; 
Objective function 2*xl+x2 is essential from step 6 

> nonessential ( [2*xl+x2 , xl+3*x2] , {xl+x2<=l , -xl-x2<=-l}) ; 
Objective function xl+3*x2 is essential from step 6 

6 Conclusion 

There are theoretical and practical reasons for developing a method to find 
nonessential objective functions. In this paper we present such a method 
and its implementation in Maple. Our algorithm is based on necessary and 
sufficient conditions for an objective function to be nonessential, and need 
only to solve a finite number of single objective linear optimization prob- 
lems. Examples showing the usefulness of our Maple package are provided: 
identification of nonessential objective functions permits to simplify the cor- 
respondent mathematical model. 
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Appendix — Maple Definitions 



Our Maple definitions follow below. The reader can download the code from 
[http : //www. mat .ua.pt/delfim/essential .html] together with many more 
examples than the ones we are able to provide in the paper. 

We begin by implementing the Gal-Leberling method (see Examples 14.11 
and 14.21) , which is based on the results of Ref . 

> ################################################################## 

> # GL method; returns true if F[-l] is nonessential; 

> # returns false if one can not conclude nothing from GL method 

> ################################################################## 

> glm := proc(F,NumVar) 

> local c, f, v, LV, lc, SolSet, SS, i; 

> c := proc(var,exp) 

> local v; 

> v:= select (has, exp+abm, var) ; 

> if v = NULL then 

> return (0) ; 

> else 

> return (v/var) : 

> f i: 

> end proc; 

> f := o -> if type (o .numeric) then o else fi: 

> v := (exp,n) -> Vector ( [seq(f (c(x I I i , exp) ) , i=l . .n)] ) : 

> LV := [seq(LinearAlgebra[VectorScalarMultiply] (v(F[i] ,NumVar) , 

alphal |i) ,i=l. .nops(F)-l)] ; 

> lc := Vector (NumVar) ; 

> for i in LV do 

> lc := LinearAlgebra [VectorAdd] (lc , i) : 

> od; 

> SolSet := solve({seq(lc [i]=c(x| I i,F[-l] ) ,i=l . .NumVar)}) ; 

> if SolSet = NULL then 

> return(f alse) ; 

> else 

> SS : = {seq(simplex [maximize] (alpha I I i , SolSet , 

NONNEGATIVE) , i=l . . nops (F) -1) } ; 

> if SS = {{}} then 

> return(f alse) ; 

> else 

> return(true) ; 

> fi: 

> fi: 

> end proc: 
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For illustrative examples on how to use the procedures stepl and step2 
see Examples EOH S3 and H2D 



> ################### 

> # step 1 

> ################### 

> stepl := proc(F) 

> local of, const, S; 

> of := add(v| 1 1 , 1=1 . .nops(F)) ; 

> const := seq(-F[i]+v| I i=0,i=l. .nops(F)) ; 

> const := [const , seq(v I I i>=0 , i=l . .nops (F) )] ; 

> S := simplex [maximize] (of , const) ; 

> return(not (evalb(subs(S ,of )=0) ) ) : 

> end proc: 

> ################### 

> # step 2 

> ################### 

> step2 : = F -> stepl (F [1 . . -2] ) : 

Follows our implementation in Maple for Step 3 (see Example 14. 6p . 

> step3 := proc(NumVar,X) 

> local LHS, RHS, SCI, SC2 , SC3, SC, SS3, i, a, epsilon, mylhs, myrhs ; 

> epsilon := 0.001; 

> mylhs := E -> if type (lhs (E) , numeric) then rhs(E) else lhs(E) fi: 

> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> LHS := [seq(mylhs(i) ,i=X)] ; 

> RHS := [seq(myrhs(i) ,i=X)] ; 

> SCI := {seq(LHS[i]+a+v| |i = RHS [i] , i=l . .nops (LHS) )} ; 

> SC2 := {seq(v||i >= , i=l .. nops (LHS) ) , a>=0>; 

> SC3 := {seq(x| | i>=0 . 001 ,i=l . .NumVar)}; 

> SC := SCI union SC2 union SC3; 

> SS3 := simplex [maximize] (a, SC) ; 

> assign(select(has,SS3,a)) ; 

> if a = then return(f alse) else return(true) fi; 

> end proc: 

In Step 4 we use a rank method, computing the rank of a matrix A by 
the Rank command from the standard LinearAlgebra package of Maple 
system. We notice that the procedure step4 does not use the last objective 
function (the last objective is given in Maple by F[-l], and we exclude it 
from consideration by writing F [1 . . -2] ). The auxiliary procedure matrixA 
is used both by Steps 4 and 7. The procedure Proposition4dotl2 is our 
Maple definition for Remark 14.41 
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> matrixA := proc(X,NumVar) 

> local c, LHS, row, mylhs ; 

> mylhs := E -> if type (lhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> c := (var,exp) -> if evalb({select(has,exp,var)} = {}) then 



else 

select (has , exp , var) /var 
fi: 

> LHS := [seq(mylhs(i)+abm,i=X)] ; 

> row := (exp, NumVar) -> map(c , [seq(x I I i , i=l . . NumVar)] ,exp) : 

> return(Matrix(map(row,LHS,NumVar))) ; 

> end proc: 
> 

> Proposition4dotl2 := proc(F,X,LE) 

> local SE, v, ETS , solW, SS, SV; 

> SE := NULL; 

> for v in LE do 

> SV := subs(v.F) ; 

> SE := SE, add(SV[i]*w| li,i=l. .nops(SV)) ; 

> od; 

> ETS := seq(SE[l]=i,i=SE[2. .-1]) , add(w I I i , i=l . .nops (SV) ) = 1 ; 

> solW := solve({ETS, add(w I I i , i=l . .nops(SV) )=1» 

union {seq(w| I i>=0 . 00001 , i=l . .nops(SV))}; 

> SS := {seq(simplex [maximize] (w| I i,solW) ,i=l. .nops(SV))}; 

> return(not(remove(i->i={},SS) = {})); 

> end proc: 

> step4 := proc (F,X, NumVar) 

> local NX , NNumVar , b , A, rnk,dif , LP, zero , efficient ,v, admissible , 

> vs, AM, sol, of ,cstEps, est, SC,p,myrhs, mylhs, LE, tv, val; 

> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> mylhs := E -> if type (lhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> tv := (x.s) -> myrhs(op(select(has,s,x))) : 

> b := X -> Vector( [seq(myrhs (i) ,i=X)] ) : 

> NX := {seq(mylhs(X[i] )+x| I (NumVar+i)=myrhs(X [i] ) , i=l . .nops(X) )}; 

> NNumVar : = NumVar+nops (X) ; 

> A := matrixA (NX, NNumVar ) ; 

> rnk : = LinearAlgebra [Rank] (A) ; 

> dif := NNumVar - rnk; 

> LP := combinat [choose] ( [seq(x I I i , i=l . .NNumVar)] , dif ) ; 

> zero := L -> {seq(i=0 , i=L) } : 

> efficient := true; 

> v := Vector ( [seq(x I I i, i=l . .NNumVar)] ) ; 

> admissible := sc -> not (member (false , {seq(evalb (myrhs (i)>=0) ,i=sc) 

> LE := NULL; 

> for p in LP while efficient do 
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> vs := subs ({seq(i=0 , i=p) } , v) ; 

> AM := LinearAlgebra[MatrixVectorMultiply] (A,vs) ; 

> sol := solve ({seq(AM[i]=b(NX) [i] ,i=l. .nops(NX))}) ; 

> if not (sol = NULL) and admissible (sol) then 

> sol := sol union zero(p): 

> of := add(epsilon I I i , i=l . .nops (F) -1) ; 

> cstEps := seq(epsilon I I i>=0 , i=l . .nops(F) -1) ; 

> est := seq(F [i] -epsilonl I i=subs(sol,F[i] ) , i=l . .nops(F)-l) ; 

> SC := {est} union {cstEps} union NX; 

> efficient := evalb( subs (simplex [maximize] (of ,SC,NONNEGATIVE) ,of)=0) ; 

> if efficient then 

> val := [seq(tv(x I I i , sol) , i=l . .NumVar)] ; 

> if not (member (false, map(i->type(i, numeric) , val))) then 

> LE := LE, {seq(x| I i=val[i] ,i=l. .NumVar)}; 

> fi; 

> fi; 

> fi: 

> od; 

> if not (efficient) then 

> return(ef f icient) ; 

> else 

> Proposition4dotl2(F[l. .-2] ,X,{LE}) ; 

> fi: 



> end proc: 

The procedure step5 makes use of an auxiliar procedure vert that receives 
three arguments: one solution given by the simplex method (denoted by 
Sol); an objective function of; and a set of constraints X. We remark that 
in Step 5 only the last objective function is used (that is given in Maple by 
F[-l] , where F is the list of all the objectives under consideration). 



> vert := proc (Sol, of ,X) 

> local v,S,LFV,LS,i,tv,Min,Max,LL,LV,gaa,ga,VFV,freeVar,Sub,LSub, 

s , delFreeVar , varSol , varOFl , varOF , Ver if ySol , aux , NX , mylhs , myrhs ; 

> mylhs := E -> if type (lhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> v := subs (Sol, of ) ; 

> varOFl := t -> op(select(i->not(type(i, numeric)) , [op(t)] )) : 

> varOF := of -> map (varOFl, {op (of )}) : 

> S := solve(of = v, varOF(of ) ) ; 

> freeVar := SS -> {seq(mylhs (i) ,i=select(i->mylhs(i)=myrhs(i) ,SS))}: 

> varSol := E -> {seq(mylhs (i) ,i=E)}: 

> LFV := freeVar(S) union (varSol(Sol) minus varOF(of)); 

> tv := (x,s) -> myrhs(op(select(has,s,x))) : 

> Min := (x,X) -> simplex [minimize] (x,X union S.NONNEGATIVE) : 

> Max := (x,X) -> simplex [maximize] (x ,X union S,NONNEGATIVE) : 
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> 



> 



> 



gaa 



LL 



LV 



[seq( [tv(x,Min(x,X)) ,tv(x,Max(x,X) )] ,x=LFV)] ; 
n -> [seq(i [j] , j=l . .n)] : 
: (n,m,L) -> if m=n then 



seq(LV(m) , i [m] =L [m] ) 



else 

seq(gaa(n,m+l , L) , i [m] =L [m] ) 
fi: 



> ga := L -> gaa(nops (L) , 1 ,L) : 

> VFV := {ga(LL)}; 

> Sub := (C1,C2) -> seq({seq(C2[j]=i[j] ,j=l. .nops(i) )},i=Cl) : 

> LSub := Sub(VFV,LFV) ; 

> delFreeVar := SS -> SS minus {seq(i , i=select (i->mylhs (i)=myrhs (i) ,SS) ) } : 

> NX : = X union {seq(i>=0,i=varSol(Sol))>: 

> VerifySol := PS -> not(member(false,{seq(evalb(subs(PS,i)) ,i=NX)})) : 

> aux := {seq(subs (s .delFreeVar (S) ) union s , s={LSub}) } : 

> return(select(VerifySol,aux)) ; 

> end proc : 



> step5 := proc(F,X) 

> local SolSM; 

> SolSM := simplex [maximize] (F[-l] ,X,NONNEGATIVE) ; 

> return(vert(SolSM,F[-l] ,X)) ; 

> end proc : 

Examples 14.111 and 14.121 illustrate the use of our Maple command step6. 

> step6 := proc(F,X) 

> local STEP5, of, cstEps, notEf f icient , sol, est, SC; 

> STEP5 := step5(F,X) ; 

> of := add(epsilon| I i,i=l. .nops(F)-l) ; 

> cstEps := seq(epsilon I I i>=0 , i=l . .nops (F) -1) ; 

> notEf f icient := true; 

> for sol in STEP5 while notEfficient do 

> est := seq(F [i] -epsilon I I i=subs (sol ,F [i] ) , i=l . .nops (F) -1) : 

> SC : = {est} union {cstEps} union X: 

> subs (simplex [maximize] (of , SC , NDNNEGATIVE) ,of) ; 

> notEfficient := evalb( subs (simplex [maximize] (of ,SC,NONNEGATIVE) ,of )<>0) ; 

> end do ; 

> return(not (notEf f icient) ) ; 

> end proc: 

Follows our Maple definition for Step 7 (see Example 14. 1 3j) . 

> step7a := proc(F,X,NumVar) 

> local b , A ,rnk, dif ,LP ,zero , efficient ,v, admissible , 

vs , AM , sol , of , cstEps , est , SC ,p , LE ,myrhs ; 



> 
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> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> b := X -> VectorC [seq(myrhs (i) ,i=X)] ) : 

> A := matrixA(X,NumVar) ; 

> rnk : = LinearAlgebra [Rank] (A) ; 

> dif := NumVar - rnk; 

> LP := combinat [choose] ( [seq(x I I i , i=l . .NumVar)] , dif ) ; 

> zero := L -> {seq(i=0 , i=L) } : 

> LE := NULL; 

> v := Vector ( [seq(x I I i , i=l . .NumVar)] ) ; 

> admissible := sc -> not (member (false , {seq(evalb (myrhs (i) >=0) ,i=sc)})) : 

> for p in LP do 

> vs := subs({seq(i=0,i=p)},v) ; 

> AM := LinearAlgebra [MatrixVectorMultiply] (A, vs) ; 

> sol := solve ({seq(AM[i]=b(X) [i] ,i=l. .nops(X))}) ; 

> if not (sol = NULL) and admissible (sol) then 

> sol := sol union zero(p): 

> of := add(epsilon| I i, i=l . .nops(F)-l) ; 

> cstEps := seq(epsilon I I i>=0 , i=l . .nops(F) -1) ; 

> est := seq(F [i] -epsilonl I i=subs(sol,F[i] ) , i=l . .nops(F)-l) ; 

> SC := {est} union {cstEps} union X; 

> efficient := evalb( subs (simplex [maximize] (of ,SC,NONNEGATIVE) ,of )=0) ; 

> if efficient then LE := LE, sol; fi: 

> fi: 

> od; 

> return ( [LE] ) ; 

> end proc: 

> # In step7b we change X 

> # Note: All inequalities must be given in the form Ax <= b 

> ############################################################## 

> step7b := proc (F,X, NumVar) 

> local TX, SEV, good, sel, mylhs, myrhs; 

> mylhs := E -> if type (lhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> TX := {seq(mylhs(X [i] )+x| I (NumVar+i)=myrhs(X [i] ) , i=l . .nops(X) )}; 

> SEV := step7a(F,TX,NumVar+nops(X)) ; 

> good := (v,nv) -> member (v, {seq(x I I i , i=l . .nv) }) : 

> sel := (es, NumVar) -> select (e->good (mylhs (e) , NumVar) , es) : 

> return (map (sel, SEV, NumV ar) ) ; 

> end proc: 
> 

> step7 := proc (F,X, NumVar) 

> local C, kernel, S7, tv, five, SD, basis, IBK, myrhs; 

> myrhs := E -> if type (rhs (E) .numeric) then rhs(E) else lhs(E) fi: 

> C := matrixA([seq(i=0,i=F[l. .-2])] , NumVar) ; 
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> kernel := LinearAlgebra [NullSpace] (C) ; 

> if kernel = {} then 

> printf ("Objective function °/ a is nonessential from step 7\n",F[-l]); 

> else 

> S7 := step7b(F,X,NumVar) ; 

> tv := (x,S) -> myrhs(op (select (has, S, x)) ) : 

> five := (Sl,S2,NumVar) -> Vector ( [seq(tv(x I I i ,S2) -tv(x I I i ,S1) , i=l . . NumVar)] ) : 

> SD := [seq(five(S7[l] ,S7[i] , NumVar) ,i=2. .nops(S7))] ; 

> basis := LinearAlgebra [Basis] (SD) ; 

> IBK := LinearAlgebra [IntersectionBasis] ( [basis, kernel] ) ; 

> if IBK = {} then 

> printf ("Objective function 7»a is nonessential from step 7\n",F[-l]); 

> else 

> printf ("X_E~7.a C X_E~°/„a from step 7\n" ,nops(F) ,nops(F)-l) ; 

> f i: 

> fi: 



> end proc : 

Our Maple procedure mm is based on the theory introduced in Ref. Q2] (mm 
stands for "Malinowska Method"). 



> mm := proc (F,X, NumVar) 

> if stepl(F) then 

> if step6(F,X) then 

> step7 (F,X, NumVar ) ; 

> else 

> printf ("Objective function %a is essential from step 6\n" ,F [-1] ) ; 

> fi; 

> else 

> if not(step2(F)) then 

> printf ("Objective function %a is nonessential from step 2\n",F[-l]); 

> else 

> if step3(NumVar,X) then 

> printf ("Objective function °/ a is essential from step 3\n",F[-l]); 

> else 

> if step4(F,X, NumVar) then 

> printf ("Objective function 7,a is nonessential from step 4\n",F[-l]); 

> else 

> printf ("Objective function °/ a is essential from step 4 

> and X_E~7 a C X_E~7a\n" ,F[-1] ,nops (F) -1 ,nops (F) ) ; 

> fi: 

> fi: 

> fi: 

> fi: 

> end proc: 



Follows the main procedure of our Maple package (see Section [5]). 
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> nonessential := proc(F,X) 

> local NumVar, y, cs; 

> cs := x -> [op(x)] [-1] : 

> y := sort(map(cs, remove(i->type(i, numeric) , 

map(i->op(i) , 

[seq(op(i) ,i=F) ,seq(op(i) ,i=X)])))) [-1] ; 

> for NumVar from 1 by 1 while not(evalb(x| I NumVar = y) 

or evalb(-x| I NumVar = y)) do od; 

> if glm(F, NumVar) then 

> printf ("Objective function °/,a is nonessential from GL method\n" ,F [-1] ) ; 

> else 

> mm (F,X, NumVar) ; 

> fi: 

> end proc: 
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